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ABSTRACT 

We present a short review of the current quasar (QSO) absorption hne constraints 
on possible variation of the fine structure constant, a = /he. Particular atten- 
tion is paid to recent optical Keck/HIRES spectra of 49 absorption systems which 
indicate a smaller a in the past |T], ^. Here we present new preliminary results 
from 128 absorption systems: Aa/a = (—0.57 ±0.10) x 10~^ over the redshift range 
0.2 < z < 3.7, in agreement with the previous results. Known potential systematic 
errors cannot explain these results. We compare them with strong 'local' constraints 
and discuss other (radio and millimetre-wave) QSO absorption line constraints on 
variations in a'^Qp and a^gpnie/mp {gp is the proton ^f-factor and rrie/mp is the elec- 
tron/proton mass ratio). Finally, we discuss future efforts to rule out or confirm the 
current 5.7 a optical detection. 
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1 Introduction 



The assumption tliat tiie constants of Nature remain constant in spacetime should 
be experimentally tested 0. Strong motivation for varying constants comes from 
modern unified theories ^, ^. Here we review the QSO absorption hne constraints 
on possible variation of the electromagnetic coupling constant, a. Our most recent 
published results are summarized in reference and in Section § we present new 
preliminary results from a significantly extended optical sample. 

2 QSO absorption systems and the alkali doublet method 

For small variations in a, the relative wavelength separation between the transitions 
of an alkali doublet (AD) is proportional to a. Savedoff 0] first utilized this to con- 
strain possible variations in a from AD separations seen in galaxy emission spectra. 
The advantage of this technique is the large look-back times inherent in such cos- 
mological observations (~ 10 Gyr). Absorption hues produced by intervening clouds 
along the line of sight to QSOs are substantially narrower than intrinsic emission 
lines and so yield tighter limits on a- variation [Q. 

Spectrographs on 8-10-m optical telescopes can record high resolution 
(FWHM ~ 7kms~^), high signal-to-noise (S/N ~ 30 per pixel) spectra of high red- 
shift QSOs over most of the optical range (i.e. 3000-8000 A) in several ~ 1 hr 
exposures. Fig. 1 shows an example QSO spectrum with a C iv AD. Many velocity 
components of the absorption system are clearly resolved. 



Varshalovich et al. |10| have recently used spectra of 16 Si IV ADs with 
a mean redshift (^abs) = 2.6 to obtain a value for the fractional difference be- 
tween a in the laboratory and in the QSO spectra, Aa/a = (a^ — Q;o)/ao = 
(—4.6 ± 4.3stat =t 1.4sys) X 10~^. The systematic error term arose from uncertainties 
in the laboratory wavelengths of the Si iv transitions: the astronomical spectra were 
of comparable quality to UV laboratory spectra. Significant improvements in the 
laboratory wavelengths []TT] reduce this systematic error to 0.2 x 10~^ (Icr)- We 



have analyzed 21 high quality Si IV doublets observed with the HIRES spectrograph 



on the Keck I 10-m telescope in Hawaii, finding [|12 | 



Aa/a = (-0.5 ± 1.3stat) X 10^^ (1) 

at (^^abs) = 2.8. The factor of 3 improvement in precision is due to the high spectral 
resolution of the HIRES data. This is currently the strongest constraint on Aa/a 
from the AD method. 
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Figure 1: Keck/ HIRES spectrum of QSO GB 1759+ 7539 j^. The full spectrum 
(upper panel) shows several emission lines intrinsic to the QSO (Ly-a, Ly-f3, NlW, 
Siw, Civ). The damped Ly-a system (DLA) at z^hs = 2.6253 gives rise to heavy 
element absorption lines in the red portion of the spectrum. The lower panel details a 
small region containing a Cw alkali doublet. The separation between corresponding 
velocity components in the two transitions is proportional to a for Aa/a <^ 1. 



3 The many-multiplet method 



The AD method is simple, but inefficient. The s ground state is most sensitive 
to changes in a (i.e. it has the largest relativistic corrections) but is common to 
both transitions (Fig. 2a). A more sensitive method is to compare transitions from 
different multiplets and/or atoms, allowing the ground states to constrain Aa/a 
(Fig. 2b). This is the many multiplet (MM) introduced in []T3|, . 

To illustrate the MM method, consider the following semi-empirical equa- 
tion for the relativistic correction. A, for a transition from the ground state with 
total angular momentum, j: 



A oc (Zaf 



3 + 1/2 



C 



(2) 



where Z is nuclear charge and many-body effects are described by C ~ 0.6. To 
obtain strong constraints on Aa/a one can (a) compare transitions of light {Z ~ 10) 
atoms/ions with those of heavy {Z ~ 30) ones and/or (b) compare s-p and d-p 
transitions of heavy elements. For the latter, the relativistic corrections will be of 
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Figure 2: (a) The AD method is not sensitive to the maximal relativistic corrections 
in the common S ground state, (h) Comparison of different ions increases sensitivity 
to Aa/a, increases statistics and decreases systematic errors. 



opposite sign which further increases sensitivity to a-variation and strengthens the 
MM method against systematic errors in the QSO spectra (see Section 

More formally, we may write the following equation for the rest-frequency, 



of any transition observed in the QSO spectra at a redshift z: 



1^0 + g 



^1 -1 

a 



(3) 



where ujq is the frequency measured in the laboratory on Earth (we omit higher order 



terms here for simplicity). Laboratory measurements |]T5|, |T6|, |TT[ of ujq for many 
transitions commonly observed in QSO spectra now allow a precision of Aa/a ~ 
10~^ to be achieved. The q coefficient contains all the relativistic corrections and 
measures the sensitivity of each transition frequency to changes in a. These have 
been calculated in [[Tsj |l7i [T^ to < 10% precision using the Dirac-Hartree-Fock 



approximation and many-body perturbation theory. Note that the form of Eq. |^ 
ensures one cannot infer a non-zero Aa/a due to errors in the q coefficients. 

Fig. 3 shows the distribution of q coefficients in (rest) wavelength space. 
Our sample conveniently divides into low- and high-z subsamples with very different 
properties. Note the simple arrangement for the low-z Mg/Fe ll systems: the Mg 
transitions are used as anchors against which the large, positive shifts in the Fe ll 
transitions can be measured. Compare this with the complex arrangement for the 
high-z systems: low-order distortions to the wavelength scale will have a varied 
and complex affect on Aa/a depending on which transitions are fitted in a given 
absorption system. In general, the complexity at high-z will yield more robust values 
of Aa/a. 
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Figure 3: Distribution of q coefficients for the transitions used in the MM method. 
For the low-z Mg/Fe systems, a compression of the spectrum can mimic Aa/a < 0. 
However, the complex arrangement at high-z indicates resistance to such systematics. 



4 Recent results 

For each absorption system we fit multiple velocity component Voigt profiles to all 
available (typically ~5) MM transitions. We minimize for all velocity components 
simultaneously to obtain the best fitting value of Aa/a. The la error is derived 
from the diagonal terms of the final parameter covariance matrix. Monte Carlo 
simulations demonstrate the reliability of both Aa/a and the errors. 

The MM method was first applied to 30 low-z Mg/Fe systems and pro- 
vided a tentative non-zero Aa/a [|14|. In we extended the sample of to 



49 absorption systems, finding 4.1 a evidence for a smaller a in the redshift range 
0.5 < 2;abs < 3.5. We have now increased our sample to 128 absorption systems, all 
observed with Keck/HIRES. Our new preliminary weighted mean is 

Aa/a = (-0.57 ±0.10) X 10"^ (4) 

for 0.2 < 2abs < 3.7, i.e. 5.7 cr statistical evidence for a smaller a in high redshift 
absorption systems. We plot Aa/a versus ^abs in Fig. 4. Note the overall internal 
consistency of the results. Fig. 5 suggests possible evolution of a with cosmological 
time, although see Section ^for further discussion and caveats of fixing Aa/a = 
at z = 0. 



5 



10 



-10 


-1 

-2 



Fractional look-back time 
0.2 0.4 0.6 0.8 



J 1 1 1 M M 1 1 

- fi 


1 1 1 1 1 1 1 

fff V\ 

1 — — \ — \ — \ — \ — — h 


H — \ — h 


1 1 _ 

It It : 

— \ — \ — \ — h 










; 1 i 

1 1 1 


1 1 1 1 1 1 


il 

1 1 1 





1 2 

Redshift 



Figure 4: Distribution of Aa/a over absorption redshift. The upper panel shows 
Aa/a for 128 absorption systems with la errors. We bin Aa/a in the lower panel, 
presenting the weighted mean Aa/a and la error at the mean redshift for each bin. 



5 Systematic errors? 

The statistical error in this result is now small: we do detect line shifts in the QSO 
spectra. But are the line shifts due to systematic errors or really due to varying a? 
We have thoroughly searched for possible systematic errors in our previous results 
finding none which provide an alternative interpretation of the data. We have 
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extended this search to the new data in Fig. 4 with similar results. Currently, our 
two largest sources of possible systematic error are: 

1. Atmospheric dispersion effects: Before 1996 Keck/HIRES had no image rota- 
tor and so the effects of atmospheric dispersion on the wavelength scale could 
not be avoided. Effective compression of the spectra may result, possibly mim- 
icking a negative Aa/a at low-^. 77 of our 128 absorption systems could have 
been affected. However, we find no evidence for these effects: the "affected" 
and "unaffected" subsamples yield the same Aa/a. Nevertheless, we modeled 
the potential effect and correct the 77 affected systems in Fig. 6 (top panel). 
This correction reduces the significance of the low- 2; points but increases the 
significance of those at high-2;, enhancing the apparent trend in Aa/a with z. 
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Figure 5: Temporal variation in a. The points are the binned values of Aa/a from 
Fig. 4> the solid line is the weighted mean and the dashed line is a fit to the raw 
(i.e. unbinned) data fixed to Aa/a = at z = 0. A analysis indicates that an 
evolving Aa /a is preferred. 



2. Isotopic ratio evolution: We fit Mg and Si absorption lines with terrestrial 
values of the isotopic ratios. If the isotopic abundances in the absorption 
clouds are different to the terrestrial values then we may introduce artificial 
line shifts, potentially leading to Aa/a ^ 0. Galactic observations and 



theoretical models strongly suggest that only the ^'^Mg and ^^Si isotopes 
will exist in the absorption clouds with significant abundances. The middle 
panel of Fig. 6 shows that the low-^ points become more significant when we 
fit only these isotopes to our data. 

The above two effects cannot explain our results. Indeed, applying both 
corrections (lower panel of Fig. 6), yields a more significant result. 



6 Other QSO absorption line methods 

Comparing absorption lines of the hydrogen hyperfine (21-cm) transition and millimetre- 
wave molecular rotational transitions offers an order of magnitude gain in precision 
over the MM method (per absorption system). The 21-cm/mm frequency ratio is 
oc o?gp Aa/a ^ manifests itself as a difference between the 21-cm and 



mm absorption redshifts. A similar difference may arise between 21-cm and optical 
absorption lines, in this case constraining a^gpfne/rnp [^. Here, gp is the proton 
gf-factor and me/rrip is the electron-proton mass ratio. 

However, systematic errors in these techniques are more difficult to quan- 
tify than for the optical MM method. Since the radio and mm continuum emission 
come from separate regions of the background QSO, the 21-cm and mm absorption 
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Figure 6: Summary of systematic effects. The top two panels compare the uncor- 
rected results in Fig. 4 with those corrected for our two most important potential 
systematic effects: atmospheric dispersion and isotopic ratio evolution. The lower 
panel combines the two corrections. 



may occur along slightly different sight-lines. Thus, a statistical sample of such 
measurements is required. Unfortunately, due to the paucity of known absorption 
systems, only two 21-cm/mm comparisons p5|, ^ and one 21-cm/optical compari- 
2^ presently exist (Fig. 7). 



son 



7 Other limits on varying a 

We summarize the strongest current constraints on a-variation in Fig. 7. For brevity, 
we do not discuss the reliability of the 'local' constraints, and refer the reader to 



27| for a review. 

Instead we focus on a comparison of the local and cosmological constraints. 
Despite the tight limits on Aa/a from laboratory atomic clocks, the Oklo phe- 
nomenon and meteoritic /3-decay, a simple non-linear evolution of a with time can 
explain all results simultaneously. Moreover, we emphasize that it is dangerous to 
compare local and cosmological limits without a better understanding of possible 
spatial variations in a 0, |2^. For example, absorption spectroscopy of ^abs ~ ab- 
sorption clouds in our Galaxy may not yield Aa/a = (cf. Fig 5). Even comparing 
the different QSO absorption constraints is difficult since the MM and AD methods 
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Figure 7: Comparison of the strongest current constraints on Aa/a from 'local' tests 
(open symbols) and QSO absorption lines (solid symbols). 



constrain a whereas the 21-cm/mm and 21-cm/optical methods constrain a'^Qp and 



a 



■gprrie/mp respectively. 



8 The future 

Ahhough the resuhs in Fig. 4 would have tremendous theoretical implications, 
confirming or refuting them is an observational issue. We are taking two main steps 
to check our recent results: 

1. Independent optical data. The greatest present concern is that only one in- 
strument has been used for all our observations. QSO spectra of similar quality 
to the Keck/HIRES data are now becoming available. Data from other tele- 
scopes/instruments (e.g. VLT/UVES) will provide an important check on our 
results. 

2. Further 21-cm/mm/optical comparisons. We are carrying out observations 
aimed at identifying new Hi 21-cm and mm-band molecular rotational ab- 
sorption systems (e.g. [^). Obtaining a statistical sample of 21-cm/mm and 
21-cm/optical comparison is vital for negating the line-of-sight velocity differ- 
ences discussed above. 

If step (1) confirms our Keck/HIRES results then step (2) will be a crucial 
check with entirely different systematic errors. 
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